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A novel effect in a granular material under swirling motion of the container is presented. At low 
packing densities the material rotates in the same direction as the swirling motion of the container 
(rotation). At higher densities the cluster of granular material rotates in opposite direction (rep- 
tation). The change of the direction of the motion of the cluster takes place at a critical packing 
density while the diffusion coefficient changes significantly. The measured critical density of the 
packing is in good agreement with results obtained by molecular dynamics simulation. 
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The fascination of the flow behavior of granular ma- 
terial like sand is due to the fact that under certain cir- 
cumstances it behaves like a fluid - an hour glass is the 
most popular example. Its solid-like behavior is obvious 
as well: A pile of sand is stable and it can undergo plastic 
deformation, but, unlike a fluid, it does not dissolve un- 
der gravity. The exotic behavior of granular matter and 
the sometimes unexpected observations have stimulated 
many scientists to focus their scientific interest to sand- 
like materials (see e.g. For a recent review see j|. 

We present a novel effect which demonstrates to a cer- 
tain amount both the similarity and the dissimilarity be- 
tween granular material and fluid. The experiment is ex- 
tremely robust, easily visualizable, and can be performed 
even as a kitchen table experiment using a bunch of mar- 
bles and a pot. The basic idea is to bring the pot in a 
swirling motion, the motion one frequently uses to stir up 
the bouquet of a glass of wine. If the pot contains only 
one or a few spheres, they will follow the rotation of the 
swirl in a similar fashion as the wine in the glass does. 
When one successively increases the number of spheres 
the angular velocity of the cluster decreases with increas- 
ing density. At a certain critical density an intriguing ef- 
fect appears: The angular velocity becomes negative, i.e. 
the cluster rotates in the opposite direction. Now the col- 
lective behavior is more reminiscent of a pancake rotated 
by the swirling motion of a frying pan. The first regime 
of motion we call "rotation" , the second one "reptation" . 
In this paper we present the first qualitative and quan- 
titative description of this crossover effect. The experi- 
ments are accompanied by a two-dimensional molecular 
dynamics simulation. The numerical results are quali- 
tatively and quantitatively in good agreement with the 
experimental findings. 

The experimental setup is shown in fig. EL A Petri dish 



of inner diameter 9.0 cm is mounted on a table perform- 
ing a circular vibration, i.e. each point of the table moves 
along a circular line during one oscillation period. For 
the experiments presented here, its frequency is chosen 
to be 2.5 s _1 and the amplitude of the swirling motion is 
1.27 cm. A camera is fixed on the oscillating table, thus 
the analysis can be done in a comoving frame by visual in- 
spection of the image displayed on a monitor. In between 
the Petri dish and the camera there is a glass plate with 
a marker which is fixed in the laboratory frame. Thus 
the position of the marker relative to the dish indicates 
the direction of the momentary acceleration. The Petri 
dish is partially filled with a monolayer of ceramic spheres 
(density l.74g/cm 3 ) of diameter 1.223 ± 0.035cm. The 
number of particles ranges between 1 and 42, where 42 
corresponds to the closest packing of spheres. We want 
to focus here only on the density dependence of the ef- 
fect. The dependence on the amplitude, the frequency, 
the material properties etc. will be discussed in more 
detail in a forthcoming paper. 

The system described above was simulated using two- 
dimensional molecular dynamics. We applied the soft 
particle Ansatz by Cundall and Strack ^ including in- 
teraction between colliding particles in normal and tan- 
gential direction. Two colliding particles i and j feel the 
force 

F l3 =F^n + F^t (1) 

with 
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Y = 8 • 10 6 g cm 2 is the Young modulus, -fN — 800 s _1 
and 7t = 3000 s" 1 are the damping coefficients in nor- 
mal and tangential direction and /i = 0.5 is the Coulomb 
friction coefficient. The unit vectors in normal and tan- 
gential direction are given by 
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FIG. 1. Experimental setup. 

The parameters are chosen to represent the material 
properties of the ceramic spheres used in the experiment 
as close as possible. Eq. (Q) describes the relative veloc- 
ity of the surfaces of the particles at the point of contact 
and eq. (||) gives the effective mass. Eq. (0) takes the 
Coulomb friction law into account, saying that two parti- 
cles slide on top of each other if the shear force overcomes 
fi times the normal force. For the integration of Newton's 
equations we applied a Gear predictor corrector method 
of fifth order (§. 

The time series of images in fig. || visualizes the mo- 
tion of 26 spheres in a counterclockwise rotating Petri 
dish at low packing density. They have been taken at 
time intervals T/5, with T being the period of the driving 
motion of the Petri dish. The evolution of the positions 
of the white spheres in figs. ||a to ||f indicates that the 
cluster of spheres also rotates in counterclockwise direc- 
tion. Figs. |J(? to show corresponding to figs. to ^/ 
snapshots from a molecular dynamics simulation using 
the same parameters as the experiment. 
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FIG. 2. Temporal behavior of a system of 26 spheres 
during one counterclockwise cycle of the swirl. The ar- 
row shows the direction of the momentary acceleration. 
The lower 6 figures show the corresponding molecular 
dynamics simulation. 




FIG. 3. Temporal behavior of a system of 37 spheres. 
Top: the experiment, bottom: molecular dynamics sim- 
ulation. 

When the number of spheres is increased the rotation 
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becomes slower. At a fairly well defined density it stops. 
Increasing the density then leads to a different kind of 
motion, which is displayed in figs. || for the case of 37 
spheres. Now the cluster rotates in the opposite direc- 
tion. This motion is now better described by the name 
reptation. Animated sequences of the motion are avail- 
able via WWW @. 

The experimental snapshots show a white arrow which 
indicates the direction of the momentary inertial force. 
It is obtained by drawing a vector from the small white 
cross through the center of the Petri dish. The small 
white cross is the marker on the glass plate which is fixed 
in the laboratory frame as mentioned above. Therefore 
it rotates in the comoving frame. This procedure shows 
that the phase shift between the crescent-shaped area 
and the inertial force is different in the rotation and the 
reptation mode. 

The most striking feature of the motion of the cluster 
of granular material is the dependence of the direction 
and the absolute value of its angular velocity on the par- 
ticle number, i.e. the packing density as presented in the 
upper part of fig. ^. To obtain the experimental data 
indicated by the solid circles, the time for a single par- 
ticle located near the wall for one complete path around 
the container is measured. The inverse of this time is 
shown for different numbers of spheres. The effect of the 
change of the sign of the rotation is thus clearly demon- 
strated. This method of measurement is not sufficiently 
well defined for a small number of particles. A particle 
located at the edge of the cluster does not remain there 
but has a tendency to migrate inwards, which makes the 
characterization of the cluster movement by tracing an 
individual particle meaningless. Therefore in the case of 
low particle numbers N < 23 the experimental method 
breaks down. 

In order to avoid those difficulties in the simulation we 
used a slightly different measurement method for the ro- 
tation velocity drawn in fig. ^ (crosses). The motion of 
the center of mass is extracted from the numerical data. 
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where (. . .}t denotes the time average, r* and v* are 
the relative position and velocity of the zth particle with 
respect to the position and velocity of the center of mass: 
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FIG. 4. The angular velocity of the cluster of parti- 
cles (top). The solid circles show the experimental data 
while the crosses show the results of a molecular dynam- 
ics simulation. For N c = 32 we observe the reversal of 
the direction of the rotation. The particle numbers ac- 
cording to figures and |^ have been marked by L and 
H, respectively. The lower part of the figure shows the 
diffusion coefficient D derived from the simulation due to 
eq. (|ll]). For N approaching N c the diffusion coefficient 
changes significantly. 

The experiment and the simulation agree quantita- 
tively, showing that the rotation velocity decreases with 
particle number. For a critical number of particles N c = 
32 we observe the reversal of the direction of the rotation. 
To obtain the experimental equivalent of eq. (||) requires 
computer aided image analysis which is currently under 
construction. In the experiment the rotation velocity is 
measured from the velocity of particles close to the wall 
of the Petri dish only, while in the simulation the velocity 
is calculated due to eq. (g) taking into account data of 
all particles. Hence the particles close to the center of 
mass having lower rotation velocity lead to lower values 
of / for the results of the simulation (s. fig. |]). One can 
get a visual impression of the described behavior in || . 

The change of the sense of rotation might be caused 
by the change of the dynamical structure of the cluster, 
which depends on the packing density. For high den- 
sities one finds that the cluster rolls like a rigid body 
along the inner wall of the Petri dish, i.e. the neighbor- 
hood relations of the particles do not change significantly 
in time. At lower density the bulk of particles behaves 
more fluid-like, i.e. the relative positions of the particles 
change quickly in time. A quantitative measure of this 
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effect is given by the diffusion coefficient D shown in the 
lower part of fig. 
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Here U denotes the surrounding of the ith particle de- 
fined by 
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In our simulation we used Rjj — 3 cm. As indicated by 
fig. ||, the diffusion coefficient changes significantly for N 
approaching N c . 

The agreement between the experimental results and 
numerical calculation is remarkable, when considering 
the fact that the theory is a two-dimensional one with 
idealized boundary conditions and particle interactions. 
Obviously, this experiment can serve as a test for the 
quality of theoretical models describing the behavior of 
granular matter, in a field where robust and reproducible 
experimental effects are not easily obtained. 
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The black and white spheres have been checked to have 
identical mechanical properties. In test series we did not 
find effects like segregation which might originate from 
different bulk or surface properties. 

The animated sequences mentioned in this paper can be 
found via 



World Wide Web under URL htt 


3://summa.physik.hu- 


berlin.de/~volkhard/swirl.html 


The following 



6 "movies" from the simulation can be accessed: a) 26 
particles recorded by a resting camera, b) camera moving 
with the dish, c) camera moving with the center of mass 
of the particles, e-g) the same for 37 particles. 
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